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Abstract 

The RGB Bayer filter system consists of a mosaic of R , G, and B filters 
on the grid of the photo sensors which typical commercial DSLR (Digital 
Single Lens Reflex) cameras and CCD cameras are equipped with. Lot of 
unique astronomical data obtained using an RGB Bayer filter system are 
available, including transient objects, e.g. supernovae, variable stars, and so¬ 
lar system bodies. The utilization of such data in scientific research requires 
that reliable photometric transformation methods are available between the 
systems. In this work, we develop a series of equations to convert the ob¬ 
served magnitudes in the RGB Bayer filter system (Rb, Gb, and Bb) into 
the Johnson-Cousins BVR filter system (Bj, Vj, and Rc )■ The new transfor¬ 
mation equations derive the calculated magnitudes in the Johnson-Cousins 
Liters ( Bjcai , V/cai, and Rccai) as functions of RGB magnitudes and colors. 
The mean differences between the transformed magnitudes and original mag- 
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nitudes, i.e. the residuals, are A (Bj — B Jca i) = 0.064 mag, A (Vj — Vj ca i) 
= 0.041 mag, and A (R c — Rccai) — 0.039 mag. The calculated Johnson- 
Cousins magnitudes from the transformation equations show a good linear 
correlation with the observed Johnson-Cousins magnitudes. 

Keywords: Data Analysis; Photometry Transformation; Johnson-Cousins 
Filter System; Bayer Filter System; Open Cluster 


1. Introduction 

Photometric systems play an important role in the quantitative photomet¬ 
ric and spectroscopic study of stars and stellar systems. The Johnson-Cousins 
UBVRI system is the most popularly used broad-band photometric system, 
and most of the existing optical data have been observed with this system 
(see Bessell, 2005). 

Lately, the chance of amateur astronomers taking parts in scientific re¬ 
search is increasing especially in fields of transient objects and variable stars. 
A large amount of the photometric data for either supernovae or variable 
stars are available with amateur CCD cameras, in which the RGB Bayer 
filter system is commonly adopted. Since most of the previously compiled 
photometric data in the astronomical community have been taken with the 
Johnson-Cousins UBVRI photometric system (Rodgers et ah, 2006; Stencel, 
2013), reliable transformation equations are necessary in order to utilize the 
data from the RGB Liter system. 

There have been many works on the transformation of RGB photometry to 
Johnson-Cousins photometry in the literature (Hoot, 2007; Loughney, 2010; 
Kloppenborg et ah, 2012; Vitek and Blazek, 2012). Especially, it was shown 
that the Liter response function of Bayer Green Liter is comparable to that 
of Johnson V Liter (see Loughney, 2010; Kloppenborg et ah, 2012). Vitek 
and Blazek (2012) transformed Bayer Liter system to Johnson-Cousins Liter 
system using color indices of stars. However, there is no equation to convert 
RGB photometry to Johnson-Cousins UBVRI photometry which reLects 
spectral types, luminosity classes and metallicities of different astronomical 
sources. 

A well-known example of the transformation between the different Liter 
systems is the photometric transformation from the Sloan Digital Sky Survey 
( SDSS ) Liter system (Fukugita et ah, 1996) to the Johnson-Cousins Liter 
system. SDSS provides homogeneous ugriz photometry for objects in a 
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large fraction of the northern sky (Karaali et ah, 2005; Rodgers et ah, 2006; 
Chonis and Gaskell, 2008; Bilir et ah, 2008, 2011, 2012; Yaz et ah, 2010). 
Fukugita et ah (1996) and Smith et ah (2002) derived an g magnitude estima¬ 
tion for Sloan photometry of the standard stars from Landolt (1992). Bilir 
et ah (2005) obtained a more accurate g magnitude transformation which 
covers both ( g — r) and (r — i ) for late-type dwarf stars. They mainly used 
color-color correlations to derive the equations. Transformation equations 
are generally dependent on metallicities and luminosity classes (Karaali et 
ah, 2005; Rodgers et ah, 2006; Bilir et ah, 2008, 2011, 2012; Yaz et ah, 2010; 
Karaali and Yaz Gokge, 2013). Recently, Ak et ah (2014) derived transfor¬ 
mation equations covering a large range of metallicities, e.g. -4 < [Fe/H] 
< 0.5 dex for giant stars. Photometric data with the DSLR RGB filter sys¬ 
tem may also be transformed, using the same principle as in the SDSS to 
Johnson-Cousins magnitude conversion. 

In this paper, we present our observations of the open clusters IC4665 and 
M52, and derive the transformation equations between the RGB filter system 
and the Johnson-Cousins BVR filter system. Circumstances of the observa¬ 
tions and the preprocessing methods are described in § 2. § 3 describes the 
process of Point Spread Function (PSF) photometry to derive final transfor¬ 
mation equations for M52. § 4 shows results of the transformation including 
the residuals, and confirms the correlations between the observed data and 
the transformed data using color-magnitude diagrams. We then apply the 
derived transformation to IC4665 to check the validity of our equations. 

2. Observations and data reduction 

2.1. Observations of M52 

M52 is a 60 ± 10 Myr old open star cluster that has a spectral type of B2 
to G2. It has interstellar reddening, E (B — V) — 0.65 mag (Bonatto and 
Bica, 2006), and has similar solar metallicity which is [Fe/H] = 0 (Tiffany, 
2008). BVR data of M52 were taken in 2014 June 27 using a 4k x 4k CCD 
camera with 15 /iin pixels with the D = 1000 mm f/7.5 RC-telescope at the 
Lemmonsan Optical Astronomy Observatory (LOAO), USA. The weather 
was clear, and the atmosphere is stable with a flux variation is only under 
2%, and the seeing disk size of ~3 pixels (~2.5"). We obtained 18 images 
in each band, with 15 sec exposures in B , 9 sec exposures in V, and 8.5 sec 
exposures in R. The field of view for the BVR data of M52 is 28' x 28'. 
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KigUXG li A linearity plot for the EOS 550D camera. The dashed line represents linearity and it fits 
the curve well up to 13000 ADU. The full-well capacity is 15300 ADU. The data were obtained by taking 
flat images. 


These data were calibrated using WEBDA 1 database (Pesch, 1960; Haug, 
1970; Wramdemark, 1976; Viskum et ah, 1997; Choi et ah, 1999; Stetson, 
2000; Pandey et ah, 2001; Maciejewski and Niedzielski, 2007). The CCD was 
cooled at -110°C. 

RGB data of M52 were taken in 2014 October 09 with an IR-cut-klter- 
removed EOS 550D (Canon Inc. 2 ) camera on the D = 300 mm f/8 RC- 
telescope at Hwasangdae Observatory (IISDO). Hongcheon-gun, Gangwon- 
do, Korea. The telescope is equipped with a field corrector, which extends 
the focal length slightly. The effective focal length of the telescope is 2600 
mm. The EOS 550D camera we have used has an APS-C CMOS sensor 
(sized 22.3 x 14.9 mm) with 4.3 /jin pixels. We selected a non-process mode 
(RAW images), and a low amplification mode (ISO = 400). The full-well 
capacity of the EOS 550D’s CMOS sensor is about 15300 ADU, while the 
linearity is guaranteed up to 13000 ADU at ISO = 400 (Figure 1). Since the 
background noise increases dramatically as battery charge decreases or after 
changing the batteries (Henden et ah, 2014), we used an AC adapter instead 
of batteries. The field of view for the RGB data of M52 is 23' x 30'. The 


1 http://www.univie.ac.at/webda/navigation.html 

2 http: //www. canon, com 
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data was obtained by taking 21 images, 120 sec exposures each at 4.7 pixels 
(~1.6") seeing. Note that the Full Width at Half Maximum (FWHM) of the 
DSLR images should be larger than 2 pixels to cover all three (Red, Green, 
and Blue) pixels with the Bayer filter set. 

2.2. Observations of 10^665 

IC4665 is a 27.7 ±1.1 Myr old open star cluster with metallicity of [Fe/H] 
= —0.03 ± 0.04 (Manzi et ah, 2008; Jeffries, 2009; Smith et ah, 2011; Lodieu 
et ah, 2011). It has interstellar reddening, E (B — V) = 0.18 mag (Prosser, 
1993). BVR magnitudes of stars in IC4665 are obtained from the WEBDA 
database (Johnson, 1954; Hogg and Kron, 1955; Alcaino, 1965; McCarthy 
and O’Sullivan, 1969; Eggen, 1971; Sanders and van Altena, 1972; Muzzio, 
1973; Neckel, 1974; Landolt, 1983a; Landolt, 1983b; Menzies et ah, 1991; 
Prosser, 1993; Menzies and Marang, 1996; de Wit et ah, 2006). 

RGB data of IC4665 were obtained in 2014 April 14. The observations 
were made with the EOS 550D camera attached to a William Optics Zenith- 
Star 80 ED APO telescope (D = 80mm, f/6.8, William Optics Corporation 3 ) 
at the Kyung Hee Astronomical Observatory ( KHAO ), Gyeonggi-do, Korea. 
The camera settings were identical to those of the M52 observations. The 
field of view ( FOV) of the system is 94' x 141'. We took 30 frames with 
15 sec exposures during a photometric night with the seeing disk size of ~5 
pixels (~8.1"). 

The details of the observations, including filters, target coordinates, field 
of views, and observing conditions are summarized in Table 1. 


Table 1: Observation log 


Date 

Observatory and Telescope 

Filters 6 

Target 

R.A. 

Dec.(J2000) 

FoV a 

Seeing (") 

2014 Apr 14 

KHAO b , 0.08m telescope 

RGB 

IC4665 

17h46ml8s 

+05°43'00" 

94' x 141' 

8.1 

2014 Jun 27 

LOAO c , lm telescope 

BVR 

M52 

23h24m48s 

+61°35'36" 

28' x 28' 

2.5 

2014 Oct 09 

HSDO d , 0.3m telescope 

RGB 

M52 

23h24m48s 

+61°35'36" 

23' x 30' 

1.6 


a Field of View 

b Kyung Hee Astronomical Observatory, http://khao.khu.ac.kr 
c Lemmonsan Optical Astronomy Observatory, http://loao.kasi.re.kr 

d Hwasangdae Observatory, 915-2 Hwasangdae-ri, Naechon-myeon, Hongcheon-gun, Gangwon- 
do, Korea 

e RGB is the Bayer RGB filter system for the Canon EOS 550D, and BVR is the 
Johnson-Cousins filter system 


3 http: //www.williamoptics.com 
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Figure 2i Throughput of Bayer R, G, and B filter for the EOS 550D (Cerny, 2013). The vertical axis 
is relative sensitivity for these filters. 


2.3. Data reduction 

The RGB Bayer filter system consists of 2 x 2 pixels mosaic RGB filter 
set that covers the sensor with Red , Green , Green, and Blue colors (Tweet 
et ah, 2009), thus we need to separate the single image into images with 
different filters. We wrote a python-based program to extract pixels that 
correspond to each color to produce R, G, and B filter images. Especially 
for Green filters, we combined two Green filter images by averaging for 
each set of RGGB filters. R, G, and B filters have center wavelength, A e ff 
= 640 nm, 530 nm, and 470 nm, respectively. Figure 2 shows throughput 
of R, G, and B filters for the EOS 550D (Cerny, 2013). Since the Bayer 
filters are directly attached on the CMOS substrate, the throughput shows 
the convolution of the filter transmissions and the quantum efficiency of the 
detector. Image combining and basic preprocessing were performed with the 
Image Reduction and Analysis Facility ( IRAF ) 4 . 

Even though the CMOS of the DSLR camera was not cooled, the average 
dark current for sufficiently short exposures (~120 seconds) was low, i.e. 1 - 
3 ADU. As a result of our dark current test, we expect that the Canon image 
processor, DIGIC 4, automatically correct dark currents. Therefore we did 


4 IRAF is distributed by the National Optical Astronomy Observatory, which is oper¬ 
ated by the Association of Universities for Research in Astronomy, Inc., under cooperative 
agreement with the National Science Foundation. 
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not apply dark subtraction to the object images. We performed bias and 
flat corrections before we break images into sub-frames, then we combined 
30 images for each filter to make a single object image to reduce the random 
noise. 

3. Results 

3.1. Photometric data 

After creating the combined images, we derived magnitudes for each star 
in the star cluster in each filter with the IRAF/DAOPHOT package. Either 
moffatl5 function or penny2 function was chosen as the PSF. We detected 564 
stars at 4cr over sky background in the M52 RGB image. We removed field 
stars that do not belong to the cluster using color-magnitude and color-color 
correlations to get only star cluster members (Pesch, 1960). After identify¬ 
ing the same stars in BVR images, we compared mag (Blue) with mag(R), 
mag(Green) with mag(V), and mag (Red) with mag(R) for individual stars to 
derive the conversion formula that describes the correlation between the two 
filter systems. Based on the color-magnitude diagram, we separate giants, 
so that only main-sequence stars are selected for transformation. All magni¬ 
tudes are de-reddened applying interstellar reddening, Ay = 3.1 xE(B — V). 
Extinction magnitudes for the RGB Bayer filter system are derived by the 
extinction law (Cardelli et ah, 1989; O’Donnell, 1994; Schlegel et ah, 1998). 
The radius of the aperture for photometry was 9.4 pixels (~3.2 W ) for the 
RGB , and 6 pixels (~5") for the BVR data. Stars in M52 have magnitudes 
and colors in the range of 8.49 < V < 14.14 mag and -0.15 < B — V < 0.48 
mag, respectively. 

3.2. Transformation equations 

The transformation equations are described as follows: 

Bj = Bb,zp + Bb + Cb,bg(Bb — Gb) + Cb,gr(Gb — Rb ) (1) 

V.j = Gb,zp + Gb + Cv : bg(Bb — Gb) + Cy^piGs — Rb) (2) 

B-c = Rb,zp + Bb + Cb,bg(B b — Gb) + Cb,gr(Gb — Bb) (3) 

In these equations, Bj, Vj, and Rc are the Johnson B , V, and Cousins R 

magnitudes, respectively. Also, Bb, Gb, and Rb indicate Bayer B, G, and R 
magnitudes. In addition, we include B b ,zp, Gb,zp, and Rb,zp representing 
the zero-points for the Bayer filters. Zero-points vary depending on the 
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observing conditions such as weather, air mass, and the system sensitivity. 
The coefficients were derived by the iterative least square fitting algorithm 
with sigma clipping method (see Table 2). 


Table 2: Coefficients and RMS variation by Sigma Clipping method 



Original Data 

1 st SC a 

2 nd SC a 

Bb,zp 

-0.371 

-0.286 

—0 ?Q1 + U - U2ei 

U.Z J1_q 0 28 

Cb.bg 

0.246 

0.419 

0-2801S 

Cb.gr. 

0.213 

0.543 

0.600d°;}?3 

RMS 

0.220 

0.083 

0.064 

Gb,zp 

-0.322 

-0.263 

0 9^0+0.022 
U.ZOZ_ 0 0 25 

Cv.bg 

0.578 

0.713 

0.542i«;i| 

Cv,GR 

-0.471 

-0.196 

-0.064i° : ^ 

RMS 

0.185 

0.057 

0.041 

Rb,ZP 

-0.277 

-0.240 

—0 99fi+ a(m 
U.ZZU_ 0 031 

Cr.bg 

0.123 

0.226 

0.0511g;|| 

Cr.gr 

0.068 

0.298 

0.468±8;it§ 

RMS 

0.163 

0.055 

0.039 


a Sigma Clipping method with 2.5a 


Figure 3 shows the distribution of magnitude residuals for the Johnson- 
Cousins systems, i.e. the difference between the observed magnitudes and the 
calculated magnitudes based on the transformation equations. Most of the 
objects fall near the zero residual line, which confirms the reliability of our 
transformation equations. However, there are a number of outliers especially 
in B band. 

4. Discussion 

4-1. Accuracy of the RGB transformation 

We derived Root Mean Square (RMS) values to test the accuracy of the 
transformation equations. The RMS errors for M52 are 0.064, 0.041, and 
0.039 mag for Bj, Vj and Rc, respectively (Table 2). 

As Figure 3 shows the dispersion around the zero line is the smallest 
for Rc while it is the largest for Bj. The RMS values also suggest that 
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Figure 3i Residuals plots for the Bj, Vj, and Rc bands for M52 from top to bottom. The horizon¬ 
tal axis indicates the observed Johnson-Cousins magnitudes. Dashed lines represent the position of lcr 
dispersions, ±0.064, ±0.041, and ±0.039 for the Bj, Vj, and Rc bands. 
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Figure 4i Diagram of RGB colors versus BVRq colors for stars in M52. The dashed lines represent 
1 st order linear fitting lines. 


the conversion into Rq magnitudes is more reliable than that of the Bj 
magnitudes. 

The distributions of the data in two color-color planes corresponding to 
two systems in Figure 4 are linear, and the dispersions are relatively small 
which indicate that our transformations are reliable. The RMS errors may 
arise from the contamination of absorption lines and from the interference 
among wavelength coverage of Bayer R , G, and B bands. 

In Figure 5, we plot the color-magnitude diagrams of M52. All of the 
stars in M52 fall onto the main-sequence. The similarity of the color magni¬ 
tude diagrams constructed using the observed Johnson-Cousins filter system 
magnitudes and the transformed magnitudes shows that the transformation 
equations work well. 


v • ; .. 



Figure 5l M52 color-magnitude diagrams from observed data (top) and calculated data (bottom). 
The dashed line represents the main-sequence (Cox, 2000). The color-magnitude diagram based on the 
magnitudes from our transformation (bottom) is consistent with that of the observed magnitudes (top). 
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Figure 6: Isochrones of M52 and IC4665 (Cassisi et al., 2006). The dashed line and cross symbols 


indicate the isochrone and observed data points of M52. And the solid line and triangle symbols indicate 
the isochrone and observed data points of IC4665. M52 and IC4665 have similar metallicities with different 
ages. 


Coefficients in our results have uncertainties derived from the reduced 
chi-squared method. 




where n is the number of stars, O is the observed data and E is the calculated 
data, and a is the error the observed data and the calculated data. Derived 
uncertainties by the Equation (4) are indicated in Table 2. 

4-2. Application to IC4665 

We applied our transformation to the RGB magnitudes of IC4665. IC4665 
is a younger open cluster, which consists of main-sequence stars and sub¬ 
giants. We excluded sub-giants for transforming stars as our transformation 
is focused on main-sequence stars. 
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Figure 7'. Residuals plots for the Bj, Vj, and Rq bands for IC4665 from top to bottom. The horizontal 
axis is the measured Johnson-Cousins magnitude. Dashed lines represent the position of la dispersions, 
±0.095, ±0.094, and ±0.070 for the Bj , Vj, and Rq bands. 
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Figure 8i IC4665 color-magnitude diagrams from the observed data (top) and from the calculated 
magnitudes (bottom). The dashed line represents the main-sequence. The calculated data were derived 
using the coefficients estimated from M52. Only the zero-points adjusted for IC4665. The color-magnitude 
diagram based on calculated data (bottom) are consistent with that of the observed Johnson B and V 
magnitudes (top). 


Figure 6 shows isochrones of M52 and IC4665. Both have similar metallic- 
ities with different ages. Since the isochrones of M52 and IC4665 are similar, 
the filter transformation coefficients can be same. The RGB magnitudes of 
stars in IC4665 were also measured using IRAF/DAOPHOT with the same 
method we used for M52. 

45 stars were detected with 4cr in main-sequence stars, and their RGB 
magnitudes were transformed into BVR magnitudes using the equations de¬ 
rived in § 3.2. All coefficients in Equations (1) - (3) are kept the same except 
the zero points, because the observing conditions of IC4665 are different from 
those of M52. We select the offset zero points at which the average residual 
becomes zero. The zero points of IC4665 are -5.813, -5.400, and -5.086 for 
B , V, and R bands. 

The RMS errors for 1C4665 are 0.095, 0.094, and 0.070 mag for Bj, Vj and 
Rc , respectively. Figure 7 shows the residuals of IC4665 between observed 
magnitudes and transformed magnitudes. 

Figure 8 shows the color magnitude diagrams of IC4665 constructed us¬ 
ing the observed magnitudes from WEBDA (top), and the Johnson-Cousins 
magnitudes calculated from the RGB magnitudes (bottom). The color- 
magnitude diagrams from the transformed magnitudes follow the same trends 
as the diagrams from the observed BVR values. The similarity shows that 
the coefficients work well not only for M52 but also for other clusters. 
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5. Summary 

Equations (1), (2), and (3) convert the magnitudes RGB Bayer filter 
system to the ones in Johnson-Cousins Elter system. The color terms, ( Bb 
— Gb) and (Gb ~ Rb) are included in the transformation equations. The 
coefficients are listed in Table 2. 

There have been a number of previous attempts for filter transformation 
(e.g. Rodgers et ah, 2006; Chonis and Gaskell, 2008; Smith et al., 2002; 
Karaali and Yaz Gokge, 2013; Ak et ah, 2014). The best instrumental fits 
before our research was obtained by Hoot (2007), who presented RMS errors 
of 0.347, 0.134 and 0.236 mag for B , V, and R bands, respectively. Our 
results are more reliable than those of Hoot (2007), because of the use of two 
color terms in the transformation equations. 

Note that our work is based on a single star cluster, which contains rel¬ 
atively young stars. Although the equations can be applied to other open 
cluster data (e.g. IC4665), the magnitude conversion would be deviated in 
other spectral types or metallicities. Further works are needed in order to 
construct more accurate magnitude transformation equations between stars, 
galaxies, supernovae, and other objects. 
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